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Active site titration by a reversible tight-binding in- 
hibitor normally depends on prior knowledge of the 
inhibition constant. Conversely, the determination of 
tight-binding inhibition constants normally requires 
prior knowledge of the active enzyme concentration. 
Often, neither of these quantities is known with suffi- 
cient accuracy. This paper describes experimental 
conditions under which both the enzyme active site 
concentration and the tight-binding inhibition con- 
stant can be determined simultaneously from a single 
dose-response curve. Representative experimental 
data are shown for the inhibition of human kallikrein. 
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The concentration of enzyme active sites IE] 0 is usu- 
ally determined in a titration experiment, using an 
irreversible or a reversible tight-binding inhibitor with 
a known apparent inhibition constant (1). Conversely, 
tight-binding inhibition constants K?* are usually de- 
termined from initial velocity measurements using a 
known concentration of the enzyme (2). Thus, the pre- 
cise measurement of either quantity. IE}* or K- P? t de- 
pends on the prior knowledge of the other. 

In this paper we show under what conditions it is 
possible to measure IE],, and K* w simultaneously, from 
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a single dose-response curve. We also describe an op- 
timal method for the determination of tight-binding 
inhibition constants in high-throughput, automated 
measurements. The method consists of fitting certain 
dose-response curves twice. In an initial least-squares 
fit the presumed enzyme concentration is held con- 
stant If the apparent inhibition constant so deter- 
mined falls into a certain range, the analysis is re- 
peated while the enzyme concentration is treated as an 
adjustable parameter. 

The merits of the method are exemplified on a 
series of Monte-Carlo simulations. We found that 
under some conditions even a relatively small bias in 
the presumed enzyme concentration {e.g., 50%) 
causes a large error (e.g., more than an order of 
magnitude) in the apparent inhibition constant. This 
result shows the need for our two-stage regres- 
sion method during automated high-throughput 
analysis of pharmaceutically important tight-bind- 
ing inhibitors. 

As an example, we show a set of laboratory data for 
the inhibition of human kallikrein. The raw data, and 
a sample computer code to perform the calculations 
proposed here, can be obtained electronically via the 
Internet at the address given in footnote 1. 

METHODS 

Materials 

Human plasma kallikrein (Athens Research and 
Technology ) and kallikrein substrate. H-Pro-Phe-Arg- 
AMC (Bachem Bioscience) were purchased from the 
indicated commercial sources. The compound RR-101. 
a competitive reversible kallikrein inhibitor, was syn- 
thesized at Axys Pharmaceuticals. 
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Experimental 

Inhibitor potency and enzyme activity measure- 
ments were performed at room temperature using a 
fluorescence microplate reader (Molecular Devices). 
Plasma kallikrein (nominal concentration 4.0 nM > was 
incubated with inhibitor RR-101 (final concentration 1, 
2. 3. 4. 6. 8. 16. and 32 nM) in 50 mM Tris IpH 7.4) 150 
raM NaCL 0.05* Tween 20, 500 fiM EDTA, and 10* 
DMSO for 30 min. Control reactions in the absence of 
inhibitor were performed in replicates of 8. The reac- 
tions were initiated by the addition of substrate (350 
uM H-Pro-Phe-Arg-AMC) and the rate of substrate 
hydrolysis was measured by monitoring the increase in 
fluorescence (excitation 355 nm, emission 460 nm) over 
a 5-min period. The commercial vendor provided an 
initial estimate of enzyme active site concentration; the 
final enzyme titer and inhibition constants were deter- 
mined by the methods described in this paper. 

Computational: Monte-Carlo Simulations 

Dose-response curves for tight-binding inhibitors 
were simulated by using the Morrison equation [1], 
where \E] Q is the active enzyme concentration, K?> is 
the apparent inhibition constant, v 0 is the control ve- 
locity observed in the absence of inhibitors <[/] 0 = 0 ) 
and v is the initial reaction velocity observed at inhib- 
itor concentration [J] 0 . 



v = v 



+ \U^Jo - [/] 0 - KD 2 + 4[E]oK?"> 
Wo • 



[1] 



Each simulated dose-response curve contained nine 
points corresponding to the following inhibitor concen- 

^n 31 ^ 2 * 50 ' 0 - 6250 > 0 1563 > 0 <>39l, 0-0099, 

0 0010, and 0.0005 uM, and a data point simulated in 
the absence of inhibitors <[/) 0 = 0). The reaction ve- 
locity in the absence of inhibitor was arbitrarily set to 
unity, Uq = 1.0. 

At each chosen value of or IE] 0 , 10,000 data sets 
were generated by superimposing normally distributed 
pseudo-random noise, with constant variance and 
standard deviation set to 5% of the maximum simu- 
lated reaction velocity. Each set of simulated data was 
subjected to nonlinear least-squares regression by us- 
ing the Levenberg-Marquardt algorithm (3). 

RESULTS 

Propagation of Systematic Errors 

In the first series of Monte-Carlo simulations we 
investigated how systematic errors in enzyme concen- 
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FIG. 1. 



Propagation of systematic errors from the active enzyme 
concentration into the apparent inhibition constant. Circles and 
error bars represent averages and standard deviations, respectively 
for each set of 10,000 fitted values of K*f \ The true value of the 
inhimtion constant was KT" = 0 01 nM. The nominal enzvme con- 
centration - 1.0 nM in all cases) was held as a 'constant 
parameter in the least-squares fit of dose-response curves to Eq [1] 
The true enzyme concentration [E] JT- was varied from 0. 1 to 2.0 nM ' 
bach dose-response curve contained nine data points, at inhibitor 
concentrations described in the experimental section. 



tration propagate into systematic errors in the appar- 
ent inhibition constant. For each simulated value of 
the inhibition constant, Kl™, 200,000 artificial dose- 
response curves were generated and subsequently sub- 
jected to nonlinear least-squares fit. 

All simulated dose-response curves contained nine 
data points, representing initial velocities measured at 
different inhibitor concentrations (see experimental 
section for the corresponding values of [7J 0 ). 

Each dose-response curve was generated by using a 
fixed Kf" for the inhibition constant (for example, 0. 1 
nM) and the correct enzyme concentration [E] J""*, in a 
subsequent least-squares fit of the same curve to Eq. 
[11, the enzyme concentration was held constant at a 
IruL^ deviated fr° m the correct concentration, 
UfrJ„ (for example, 1.0 nM instead of 0.9 nM). Be- 
cause of this systematic error in enzyme concentration 
the fitted inhibition constant K\* x) was different from 
the simulated inhibition constant Kl^. 

At each combination of [E] 0 and JTf», 10,000 slightly 
different dose-response curves were simulated and 
subsequently fit to Eq. [1], yielding a spread of values 
Iu r tte ? inhibition constant. The distribution of 
the K, values was approximately Gaussian. There- 
fore, it was possible to summarize the results of 10 000 
simulations in the form of an average and a standard 
deviation. One series of numerical experiments, ob- 
tained for Kr* = 0.1 nM, is summarized in Fig. 1 
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FIG. 2. Uncertainty in the determination of K?* from dose- 
response data. Dose-response curves were simulated at [EYF"*' 
1.0 nM in all cases. For each panel a through d, 10,000 slightly 
different dose-response curves were simulated and subsequently 
subjected to least-squares fit to Eq. [1). The nominal enzyme 
concentration was treated as an adjustable parameter. Note the 
non-Gaussian character of the histogram in panel d; approxi- 
mately 20% of all data sets in that series lead to the value of K\ tu 
which was practically zero. 



Figure 1 represents 20 sets of simulations assuming 
the true value K 1 *"* = 0.1 nM and the nominal enzyme 
concentration t^Jo 000 " *= 1.0 nM, which was held con- 
stant in the least-squares fit. The true enzyme concen- 
tration [E] i,^ was varied between 0.1 and 2.0 nM. The 
results show that comparatively small systematic er- 
rors in enzyme concentrations propagate into large 
systematic errors in inhibition constants. For example, 
if the true enzyme concentration were 0.5 nM instead 
of the assumed 1.0 nM, the inhibition constant would 
be underestimated on the average by a factor of 9 (K'^ 
= 0.012). Even a systematic error as small as 10% in 
enzyme concentration ( [E] J™" - 0.9 nM) leads to fitted 
values of inhibition constants that are underestimated 
on the average by a factor of 2 to 3. 

Even more prominent distortion in the fitted value of 
K?* was seen in a similar series of experiments at 
= 0.01 nM. under otherwise identical conditions 
([E]^ 00 " = 1.0 nM). In that series of Monte-Carlo sim- 
ulations (data not shown), at some values of the true 
enzyme concentration (e.g., at [Eli'™*' = 0.1 nM) a 
significant number of the dose-response curves (ap- 
proximately 2000 of the 10.000 total data sets) could 



not be fit to Eq. [11 at all. In those cases the best-fit 
value K I*™ 4 ' was practically zero. This is illustrated also 
in Fig. 2d below. 

Simultaneous Determination ofK? p and [E] 0 

In the second series of Monte-Carlo simulations, we 
examined the conditions under which it is possible to 
determine both the inhibition constant and the enzyme 
active site concentration from a single dose-response 
curve. This simultaneous determination is based on 
fitting dose-response data to Eq. [11, while treating not 
only and v 0 but also [E) 0 as adjustable parame- 
ters. 

The simulated dose-response curves again consisted 
of nine data points (zero inhibitor concentration plus 
eight data points at inhibitor concentrations described 
in the experimental section). The simulated enzyme 
concentration was always 1.0 nM, while the inhibition 
constant was varied from 10 nM to 0.01 nM, stepping 
by a factor of 10. At each value of the apparent inhibi- 
tion constant, 10,000 slightly different dose-response 
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FIG. 3. Uncertainty in the determination of from dose-re- 
sponse data. For each panel a through d, 10,000 slightly different 
dose-response curves were simulated and subsequently subjected to 
least-squares fit to Eq. (1]. The nominal enzyme concentration was 
created as an adjustable parameter. Note the non-Gaussian charac- 
ter of the histogram in pane] a; most data sets in that series lead to 
the value uf \E]{" which was practically zero, which means that the 
enzyme concentration cannot be estimated under these conditions. 
However, the inhibition constant is estimated verv reliablv (see 
Fig. 2aK 
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FIG. 4. Inhibition of human kallikrein. The initial velocities 
(relative fluorescence units per second) were fit to Eq. [1J. Dashed 
curve: enzyme concentration was kept constant at the nominal 
value [EU = 4.00 nM; the fitted value of the inhibition constant 
was A - 8o pM. Solid curve: enzyme concentration was treated 
as an adjustable parameter; the fitted values were [El 0 = 1.53 nM 
and A w = 478 pM. 



curves were simulated. Each simulated curve was sub- 
sequently fit to Eq. [1] and the best-fit values of K?" 
and \E} 0 were recorded. The spread of best-fit values of 
K!» is shown in Fig. 2, while the spread of best-fit 
values of [E) 0 is shown in Fig. 3. 

Figure 2 shows that the histogram of distribution for 
the fitted values of K?> is reasonably narrow only for 
those values that are larger than 10% of the true en- 
zyme concentration. On the other hand, Fig. 3 shows 
that the enzyme concentration can be determined with 
certainty only if the inhibition constant is, at most, 
equal to the true enzyme concentration itself. Thus' 
between KJ» > 0.1 x [E] 0 and Kf < 1.0 x [E] 0t both 
quantities can be determined simultaneously from the 
same dose-response curve. 

Sample Experimental Data 

Human kallikrein was assayed at the nominal con- 
centration [E] 0 = 4.00 nM, based on the nominal con- 
centration of the enzyme stock. Protease inhibitor RR- 
101 was used to inhibit the enzyme in 96-well plate 
reader assays. The initial reaction velocities at differ- 
ent inhibitor concentrations are shown as open circles 
in Fig. 4. 

The Morrison equation [1] was used to fit the exper- 
imental data in two different ways. In the first round of 
analysis, represented by the dashed curve in Fig 4 the 
concentration of the enzyme was held constant at its 
nominal value [E] 0 - 4.00 nM. The best-fit value of the 
apparent inhibition constant so obtained was 85 ± 97 
pM. The relatively large standard error is caused by 



the systematic discrepancy between the experimental 
data and the theoretical model (the dashed curve in 
Fig. 4 deviates significantly from the data points). The 
sum of squared deviations was 0.3798. 

In the second round of analysis, represented by the 
solid curve in Fig. 4, the enzyme concentration was 
treated as an adjustable parameter yielding the best-fit 
value of [E] 0 = 1.53 ± 0.19 nM. This value represents 
approximately 40# of the nominal enzyme concentra- 
tion. The best-fit value of the apparent inhibition con- 
stant was 478 z 57 pM, which is approximately six 
times the value observed in the first round of nonlinear 
regression. The sum of squared deviations between the 
experimental data and the theoretical model was 
0.0141. 



DISCUSSION 

Laboratories engaged in high-throughput kinetic 
characterization of enzyme inhibitors have an increas- 
ing need to develop primary and secondary (follow-up) 
screening assays in which the concentrations of the 
enzyme and the inhibitor are comparable in magnitude 
("tight binding*). The automated two-stage regression 
analysis described in this paper represents a new data- 
analysis technique that can be used to improve the 
accuracy of such high-throughput enzyme kinetic anal- 
ysis without any additional expenditures in time or 
materials. 

We have recently described a new method for auto- 
mated determination of tight-binding enzyme inhibi- 
tion constants (4). The method has been applied suc- 
cessfully to more than 100,000 inhibitors of mast-cell 
tryptase and other therapeutic targets. In the course of 
this work, we noticed that the inhibition constants 
seemed remarkably sensitive to the presumed enzyme 
concentration. Particularly for extremely tight binding 
inhibitors, it appeared that small variations in the 
enzyme concentration produced large variations in the 
observed inhibition constant. To our knowledge, no 
systematic statistical treatment of such error propaga- 
tion exists in the analytical literature. 

In this study we present the results of Monte-Carlo 
simulations predicting that for very tight-binding in- 
hibitors (Kr - 0.1 x [E] Q or 0.01 x [E] Qt see Fig. 1), 
a systematic error in enzyme concentration as small as 
10-20% might produce as much as an order of magni- 
tude systematic distortion of the apparent inhibition 
constant. This prediction was confirmed on experimen- 
tal data for the kallikrein inhibitor RR-101. At the 
presumed enzyme concentration [E] 0 = 4.00 nM the 
mhibition constant was 85 pM, while at the best-fit 
enzyme concentration [E} 0 = 1.53 nM, the apparent 
inhibition constant was 478 pM. 
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Systematic errors in the concentration of enzyme 
active sites are probably quite common in biochemical 
laboratories, including high-throughput screening es- 
tablishments. Enzyme preparations easily change spe- 
cific activity due to a host of denaturation processes or 
small changes in the assay conditions (pH, tempera- 
ture, ionic strength). The question is how to avoid the 
consequences, in the form of large systematic errors in 
the apparent inhibition constants, of such inadvertent 
changes in total enzyme concentration. 

The strategy proposed in this report emerged from 
the results of Monte-Carlo simulations shown in Figs. 
2 and 3. We found that when the inhibition constant is 
lower than the total enzyme concentration, it is possi- 
ble to extract both quantities from a single dose-re- 
sponse curve. On the other hand, when the inhibition 
constant is relatively high (K? p > [E] 0 ), the enzyme 
concentration cannot be extracted from the experiment 
and therefore must be treated as a constant. Fortu- 
nately, under such "loose binding" conditions even a 
very large systematic error in enzyme concentration 
has no effect on the observed inhibition constant, be- 
cause the Morrison equation [11 reduces to Eq. [2]. 

Based on these results, the proper computational 
algorithm for high-throughput analysis of enzyme in- 
hibitors is the two-stage regression analysis proposed 
in this paper. In the first round of analysis, the pre- 
sumed enzyme concentration is treated as a constant 
parameter. The apparent inhibition constant so deter- 
mined is accepted as the final result if it is higher than 
the nominal enzyme concentration. Otherwise the 
least-squares fit of the same data to Eq. [11 is repeated, 
while the enzyme concentration is treated as an adjust- 
able parameter. 

Based on Monte-Carlo simulation studies similar to 
those reported here, Greco and Hakala (2) proposed 
that the accuracy of measuring tight-binding inhibition 
constants can be improved by lowering the enzyme 
concentration, if possible even to the point where tight- 
binding is avoided (that is, where [E) 0 « K* w ). How- 
ever, for many enzyme systems this strategy is not 
possible or practical, because commonly used synthetic 
substrates are too inefficient. A relatively large amount 
of enzyme would then be required to complete the 
assay in a reasonably short time. For example, assays 
of viral proteases such as cytomegalovirus assembling 
(5, 6) may require very high concentrations, in the 
micromolar range, in order to detect substantial pro- 
teolytic activity. Some mammalian proteases that are 



potential therapeutic targets such as prostate specific 
antigen (7) or urokinase (8) are routinely assayed at 
concentrations approaching 100 nM. Under such con- 
ditions even inhibitors with Kt" ** 10 nM behave as 
"tight binding." 

The two-stage regression method described in this 
paper has two important limitations. First, the nomi- 
nal enzyme concentration must not be too inaccurate. 
In exploratory Monte-Carlo simulations (data not 
shown), we found that the nominal enzyme concentra- 
tion must be accurate approximately within an order of 
magnitude. Outside this range, the nonlinear least- 
squares regression might produce an erroneous esti- 
mate of the apparent inhibition constant, ending in a 
loca minimum on the least-squares surface. 

The second limitation of the two-stage regression 
method is illustrated by Fig. 2d. A significant number 
of the simulated dose-response curves (approximately 
2000 or 20% from the total) produced artificially low 
estimates for the inhibition constants, practically ap- 
proaching zero. This problem can be minimized by a 
proper design of the kinetic experiment, that is, by 
optimally choosing the inhibitor concentrations. Our 
preliminary results suggest that an optimally designed 
series of inhibitor concentrations includes several 
points near the nominal enzyme concentration. The 
inhibitor dilution series chosen in this study (see ex- 
perimental section) is not optimal in that regard. The 
results of a detailed study on the optimal design of 
high-throughput kinetic experiments will be reported 
elsewhere. 

It has been said (9) that "in the 100 years to 1995, the 
pharmaceutical industry worked on about 500 targets 
with a limited number of compounds, whereas now, 
using new technologies like genomics, high throughput 
screening, and combinatorial chemistry, drug compa- 
nies will see an explosion in the number of targets and 
leads [they] can explore. To improve the transition 
from research to development it is necessary to auto- 
mate the research and development process [ ... 1 us- 
ing information technologies to make better use of ex- 
isting data." Our two-stage regression method 
represents a step toward this goal. A suitable computer 
program for performing the kinetic analysis is 
available free of charge via the internet (httpV/www. 
biokin.com). 
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